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An octopus program is demonstrated to generate electron energy levels in three-dimensional 
geometrical potential wells. The wells are modeled to have shapes similar to cone, pyramid and 
truncated-pyramid. To simulate the electron energy levels in quantum mechanical scheme like the 
ones in parabolic band approximation scheme, the program is run initially to find a suitable electron 
mass fraction that can produce ground-state energies in the wells as close to those in quantum dots 
as possible and further to simulate excited-state energies. The programs also produce wavefunctions 
for exploring and determining their degeneracies and vibrational normal modes. 



INTRODUCTION 



II. SYMMETRIC POTENTIALS 



Symmetry plays an important role in quantum 
mechanics^ It helps predict and classify wavefunctions 
a quantum system can have in each energy level. For 
three-dimensional rotational invariant potentials such as 
in a harmonic oscillator and a hydrogen atom, their 
Schrodinger's equations, 



Before doing quantum mechanical simulation for three- 
dimensional potential wells, it would be nice to test the 
precision of octopus with well-known symmetric poten- 
tials. Let one try first with one-dimensional harmonic 
oscillator potential, 
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can be solved analytically and there exist explicit formu- 
las which relate a number of wavefunctions or degeneracy 
to energy level n. For other symmetric, geometrical po- 
tentials, analytical method may not be able to find their 
solutions. Except by guessing from symmetry of the po- 
tentials, one hardly knows their approximate solutions. 
However, with the inclusion of computational tools, their 
numerical solutions including their degeneracies can be 
seen explicitly. 

As it is realized that simulation is one of the paths 
toward scientific truth. 3,4 In this note, programs in octo- 
pus code 5,6 are demonstrated to generate electron energy 
levels in three-dimensional geometric potential wells. Al- 
though octopus package is designed to solve many-body 
systems based on time-dependent density functional the- 
ory, it is able to solve single-particle systems based on 
the Schrodinger's equation and to produce output files 
for potentials and wavefunctions, which can be viewed 
by visualizers such as Gnuplot, OpenDX, etc. The po- 
tential wells are modeled to have shapes similar to cone, 
pyramid and truncated-pyramid^^ii^ To simulate the 
electron energy levels in quantum mechanical scheme like 
the ones in parabolic band approximation scheme, the 
programs are run initially to find ground-state energies 
in the wells as close to those in quantum dots as possible 
and further to simulate excited- state energies. The pro- 
grams also produce wavefunctions, which are visualized 
by OpenDX, for exploring and determining their degen- 
eracies and vibrational normal modes. 



and one-dimensional symmetrical linear potential, 
V(x) = \x\. 
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A bonus from the symmetric linear potential is that one 
is able to derive solutions for the linear one, 
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by symmetry consideration. Although the linear poten- 
tial is simple, one needs to look for zeros of Airy function 
to get energy eigenvalues^^ 



A. The octopus input 

For demonstration in this note, the octopus source 
code version 3.0.1— is downloaded and complied by gfor- 
tran on 64-bit Xubuntu 8.04 LTS running on Intel E7200 
with 4GB RAM. The input file for a particle called "Q- 
ball" in the harmonic oscillator potential (|3]) is shown 
below. 

(1) °/ CalculationMode 

(2) gs I unocc 

(3) % 

(4) NumberUnoccStates = 9 

(5) TheoryLevel = independent .particles 

(6) EigenSolverMaxIter = 5000 

(7) Dimensions = 1 

(8) radius = 10 

(9) spacing =0.1 
(10) °/ Species 
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(11) " Q-ball "111 user.def ined 1 1 1 " . 5*x~2 " 

(12) % 

(13) °/ Coordinates 

(14) "Q-ball" I I I 

(15) °/ 

(16) Output = wis + potential 

(17) OutputHow = axis_x + gnuplot 

(18) Output WfsNumber = "1-10" 

Line numbers are associated with the following refer- 
ences. 

• Line 1-3 form a block of calculation mode. The 
program will run initially for a ground state and 
further for the unoccupied states. 

• Line 5-6 declare type of theory and a maximum 
iteration number in solving for each eigenvalue. 

• Line 7-9 declare a dimension, a half-length of com- 
putational box and a spacing number or mesh size. 

• Line 10-12 form a block defining a particle called 
"Q-ball" having one electron mass and charge in a 
user-defined potential 0.5a; 2 . 

• Line 13-15 form a block defining Q-ball's coordi- 
nates. 

• Line 16-18 declare output options. 

For the symmetric linear potential, just replace the har- 
monic oscillator potential with abs(x). After discretizing 
the Hamiltonian over mesh points, eigenvalues are solved 
by using conjugate gradient method. Tolerance of each 
eigenvalue is of order 10 -6 by octopus' default setting. 



B. The octopus output 

Since octopus is based on real-space calculation, a 
spacing number between two nearest-neighbored mesh 
points is one of the key parameters. So, before running 
the aforementioned input, it needs to check whether the 
chosen spacing is able to produce precise energy eigen- 
values. The ground-state energies of both the harmonic 
oscillator and symmetric linear potentials varying with 
spacings are shown in Fig. [TJ Spacings starting from 
0.2 give ground- state energies of the harmonic oscilla- 
tor gradually falling off from their exact value of 0.5 
(h = uj = 1), while those of the symmetric linear po- 
tential rapidly falling off from their upper bound value. 
It is noted that for the harmonic oscillator the spacing 0.1 
produces the exact energy eigenvalues up to n — 8 and 
slightly off the exact value when n > 8. Energy eigen- 
values of both the harmonic oscillator and the symmet- 
ric linear potentials from octopus are shown in Table H 

Energy eigenvalues of the linear potential are derived 
from those of the symmetric one by imposing a bound- 
ary condition that wavefunctions must be zero at x = 0. 
Hence, only energy eigenvalues of the symmetrical linear 
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FIG. 1: Ground-state energy vs spacing of (a) harmonic os- 
cillator potential and (b) symmetric linear potential. 
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TABLE I: Energy eigenvalues of harmonic oscillator poten- 
tial, symmetrical linear potential, linear potential and zeros 
of Airy function. When the energy eigenvalues of the lin- 
ear potential are multiplied by v 7 ^, the resulting values are 
approximately equal to zeros of the Airy function. 13 



potential at even- number level are those of the linear one, 
shown in the third column of Table H 

Now, the input can be extended from one-dimensional 
to three-dimensional by setting 'Dimensions = 3' and 'ra- 
dius = 8' and replacing the potential with 

0.5*(x~2+y~2+z~2) 

Variation of the ground-state energies with spacings of 
three-dimensional harmonic oscillator has the same pat- 
tern as that of the one-dimensional; the energy eigenval- 
ues gradually fall off its exact value of 1.5 as spacings 
increase. 



III. GEOMETRICAL POTENTIAL WELLS 

Next, we will do simulation for geometrical potential 
wells. Their shapes shown in Fig. [2] are a cone of radius 
10 nm and height 10 nm^ a pryramid of width 12.4 
nm and height 6.2 nm^ and a truncated- pyramid of 
width 12.4 nm and height 3.1 nm. 9 
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FIG. 2: Geometric potential wells in the shape of (a) cone, 
(b) pyramid and (c) truncated-pyramid. These pictures are 
produced by octopus and viewed by OpenDX. 



A. The octopus input 

An input file for octopus to simulate electron energy 
levels and wavefunctions in a geometric potential well 
V(x,y,z) with depth offset of —0.77 is shown below. 

(1) Units = ev_angstrom 

(2) CalculationMode = gs 

(3) # CalculationMode = unocc 

(4) # NumberUnoccStates =11 

(5) TheoryLevel = independent_particles 

(6) # EigenSolverMaxIter = 5000 

(7) Dimensions = 3 



(8) BoxShape = parallelepiped 

(9) lsize = 150 

(10) spacing = 3.785 

(11) ParticleMass = 0.036 

(12) °/ Species 

(13) "QD" | ParticleMass I user.defined | 1 | "V(x,y,z) " 

(14) % 

(15) °/ Coordinates 

(16) "QD" | I I 

(17) °/ Discussions 

(18) Output = potential + wfs 

(19) OutputHow = dx 

(20) OutputWfsNumber = "1-12" 

The 'lsize' and 'V(v,y,z)' in the input are replaced ac- 
cording to the following potential wells: 

• Cone 

Set 'lsize' in the following block form 
°/ Lsize 

140 I 140 I 80 

7o 

and replace 'V(x,y,z)' with 

-0 . 77*step (50-abs (z) ) 
*step(50-sqrt (x~2+y~2)-z) 

• Pyramid 

Set 'lsize' in the following block form 

°/ Lsize 

124 | 124 | 60 

7o 

and replace 'V(x,y,z)' with 

-0 . 77*step(62-abs(x) ) *step(62-abs (y) ) 
*step(31-abs (z) ) *step(31-z-abs (x) ) 
*step(31-z-abs (y) ) 

• Truncated-pyramid 

Replace the pyramid potential with 

-0 . 77*step (62-abs (x) ) *step (62-abs (y) ) 
*step(31+z)*step(31-z-abs(x)) 
*step(-z) *step(31-z-abs (y) ) 

Due to regular grid discretization in all directions, there 
is only one spacing parameter. By using octopus' de- 
fault settings, eigenvalues are solved by conjugate gradi- 
ent method and their tolerance is of order 10 -6 . 



B. The octopus output 

The input files are run in the following steps: 
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FIG. 3: Ground- state energy vs spacing. 



FIG. 5: Eigenvalue distributions in the cone potential, the 
pyramid potential and the truncated-pyramid potential. 
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FIG. 4: Ground-state energy vs electron mass fraction. 



1. Explore how the ground-state energies vary with 
spacings. To run the input automatically, use a 
shell script below. 

# ! /bin/bash 
list="2 3 4 5 6 7 8 

9 10 11 12 13" 
export 0CT_PARSE_ENV=1 
for 0CT_spacing in $list 
do 

export 0CT_spacing 

octopus >& out-$0CT_spacing 

energy='grep Total static/info \ 

I head -1 I cut -d "=" -f T 
echo $0CT_spacing $energy 
rm -rf restart 
done 

Unexpectedly, when spacings increase, ground 
state energies in the wells vary in sawtooth pat- 
terns with increasing amplitudes and period widths 
as shown in Fig. [3j It seems that for each potential 
well there are many spacing parameters that can 
produce the same ground-state energy. To simulate 
the energy eigenvalues as those in Refs.iandi, the 
spacing parameter of 0.3875 nm is chosen. Due to 
time limit in class for demonstration, higher spac- 
ing parameters such as 0.775 nm is recommended. 



2. Find a suitable electron mass fraction that can 
produce ground-state energies in the geometrical 
potential wells as close to those in the quantum 
dots as possible. §a2*lfi To run them automatically 
by using the previous script, replace numbers in the 
'list' with the ones as shown on the horizontal axis 
of Fig. [4] and 'spacing' with 'ParticleMass' every- 
where. As shown in Fig. [H the ground-state ener- 
gies smoothly decrease as mass fractions increase. 
It is estimated that the suitable mass fraction is 
0.036. 

3. Explore the excited- state energies and wavefunc- 
tions by choosing the 'unocc' mode. The distri- 
butions of the energy eigenvalues in the geomet- 
rical potential wells are shown in Fig. [5l The 
wavefunctions are exported as DX-formatted two- 
dimensional cross sections for exploring their nor- 
mal modes: 

• cone 

Cross sections of wavefunctions in accordance 
with energy eigenvalues Ai, A2, A4, A5, A7, Ag 
and A10 are shown in Fig. [6l The wavefunc- 
tions of A3 and A9 are obtained from those 
of A2 and As by rotating around the i-axis 
90° clockwise and counter-clockwise, respec- 
tively, and the ones of A6 and An from those 
of A5 and A10 by rotating around the i-axis 45° 
clockwise and counter-clockwise, respectively. 
It appears that all two-dimensional cross sec- 
tions of wavefunctions in the cone look similar 
to the vibrational normal modes of a circular 
membrane 14 
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where subscript (n, m) denotes a vibrational 
normal mode of the wavefunction with har- 
monic order mth of the nth-order Bessel func- 
tion of the first kind J n (kmp)- The normal 
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FIG. 6: Cross sections of (a) ^ Ai , (b) ^a 2 , (c) ^a 4 , (d) ^a 5 , 
(e) ^a 7 , (f) ^a 8 and (g) ^A i0 in the cone potential at z = —31. 



mode (n, m) of the circular membrane is in 
general doubly degenerate except (0, m) which 
is non-degenerate. Hence, the normal mode 
of the wavefunction in the cone ^a* matches 
exactly one-to-one with that in the circu- 
lar membrane ^( n>m ). It yields the following 
equivalence in normal modes of the wavefunc- 
tions: 

*Ai ~ *(0,1)> ^A 2 , 3 ~ *(1,1), ^A 4 ~ * (0,2), 
^A 5 , 6 ~ *(2,1), ^A 7 ~ *( 0> 3), ^A 8 , 9 ~ *(l,2), 
*A 10 ,ii ~*(3,1)- (6) 



• pyramid and truncated-pyramid 
Cross sections of wavefunctions in accordance 
with energy eigenvalues Ai and A2 in the pyra- 
mid and the truncated-pyramid are shown in 
Fig.[7land[Hl respectively. Both wavefunctions 
of A3 are obtained from those of A2 by rotating 
around the i-axis 90° counter-clockwise and 
clockwise, respectively. It also appears that 
all two-dimensional cross sections of wave- 
functions in the pyramid and the truncated- 
pyramid look similar to vibrational normal 
modes of a square membrane 14 
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FIG. 7: Cross sections of (a) and (b) ^a 2 i n the pyramid 
potential at z = —15.5. 
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FIG. 8: Cross sections of (a) and (b) ^a 2 i n the 

truncated-pyramid potential at z = —15.5 . 



^(n,m) ~ sin(2ranr/a) sin(2ra7n//a). (7) 

Matching normal modes of wavefunctions in 
the pyramid and the truncated-pyramid with 
ones in the square membrane is straightfor- 
ward only for the ground state. For their 
excited- states, it appears that the nodal line of 
A2 in the pyramid is along vertical with a bit 
distortion, while one in the truncated-pyramid 
is along diagonal. These patterns are due to 
the hybridization or superposition of the fun- 
damental normal modes (1,2) and (2,1) of 
the square membrane. By configuring out a 
combination of ^(1,2) an d ^(2,1) > it yields the 
following equivalence in normal modes of the 
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wavefunctions: 
For the pyramid, 

*Ai~*(l,l), *A 2 ~-*(1 1 2)+4*( 2 ,1), 

^A 3 ~4* ( i j2 ) +*( 2> 1). (8) 

For the truncated-pyramid, 

*Ai~*(l,l), *A 2 ~-*(l,2)-*(2,l), 

^A 3 ~*(1,2) -*(2,1)- (9) 

IV. CONCLUSIONS 

Octopus is used to simulate the energy levels and wave- 
functions in the geometrical potential wells. Although 
one can guess from the symmetry of the wells what the 
possible solutions including their degeneracies would look 
like, octopus is able to verify these presumptions. Since 
octopus is based on a real-space computation, it needs 
to take a close look at the effect of the spacing parame- 
ter to the energy eigenvalue before doing simulation. For 
the cone, pyramid and truncated-pyramid, there seems to 
have many spacing parameters producing approximately 
the same eigenvalue. So, spacing of 0.3875 nm is chosen 
as in Refs. [s| and El To simulate the electron energy levels 
in quantum mechanical scheme like the ones in parabolic 



band approximation scheme, the input files are run ini- 
tially to find the electron mass fraction whose suitable 
value is 0.036. It yields the energy eigenvalues in the 
cone off from those in Ref. |l0| in an order about 0.02- 
0.05 and ones in both the pyramid and the truncated- 
pyramid off from those in in Refs. [H and in an order 
about 0.01. Octopus is also able to produce the DX- 
for matted wavefunctions. When looking at their cross 
sections, normal modes in the cone match exactly one- 
to-one with those in the circular membrane. In the pyra- 
mid and the truncated-pyramid, only normal mode in the 
ground state matches with that in the square membrane, 
while ones in their first excited-states are from the hy- 
bridization or superposition of fundamental modes ^(1,2) 
and ^(2,1) in the square membrane. Without the simu- 
lation, these vibrational normal modes might be beyond 
imagination. 
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